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ABSTRACT 


*  Many  military  communication  systems  are  required  to  operate  in  environments 

where  an  enemy  may  attempt  to  interfere  with  the  conununication.  One  of  die  most  serious 
threats  is  the  use  of  jammers.  This  threat  has  caused  designers  of  communications  systems 
to  counter  with  the  use  of  spread  spectrum  techniques  such  as  frequency-hopping. 

Frequency-hopped  spread  spectrum  satellite  communication  systems  provide  robust 
communications  resilient  to  many  types  and  levels  of  jamming.  Unfortunately  this  increase 
in  resilience  is  offset  by  an  increase  in  complexity  while  establishing  the  communication 
link,  termed  synchronization.  This  document  describes  a  downlink  synchronization  process 
that  involves  the  transmission  of  synchronization  hops  by  the  satellite  and  a  two-step 
ground  terminal  synchonization  procedure.  In  addition  computer  simulation  results  are 
documented  which  demonstrate  that  the  proposed  synchronization  algorithm  is  able  to 
operate  in  environments  as  severe  as  the  communication  link. 


R^SUMg 


Beaucoup  de  syst&mes  militaires  de  communication  sont  requis  pour  opdrer  dans 
des  milieux  oh  un  ennemi  peut  essayer  d'interfdrer  avec  les  communications.  Une  des 
menaces  les  plus  sdrieuses  est  Tutilisation  de  brouilleurs.  Cette  menace  a  forcd  la  defense 
des  systhmes  de  communication  en  utilisant  des  techniques  d'dtalement  du  spectre  telles  que 
le  saut  en  frequence. 

Les  systhmes  de  communication  par  satellite  avec  dtalement  du  spectre  par  saut  de 
frequence  foumissent  des  communications  robustes,  protdgdes  contre  plusieurs  types  et 
niveaux  de  brouillage.  Malheureusement,  cette  protection  additionnelle  est  aux  ddpens  d'une 
complexity  accrue  de  la  synchronisation,  laquelle  marque  le  ddbut  du  lien  de 
communication.  Ce  document  ddcrit  un  processus  de  synchronisation  du  lien  descendant 
impliquant  la  transmission  de  sauts  de  synchronisation  par  le  satellite  et  une  procddure  de 
synchronisation  h.  deux  dtapes  pour  la  station  terrestre.  De  plus,  des  rdsultats  ^  simulation 
par  ordinateur  sont  inclus,  lesquels  d^ontrent  que  I'algoritlune  proposd  de  synchronisation 
est  capable  d'opdrer  dans  des  conditions  similaires  k  celles  anticipdes  pour  le  lien  de 
communication. 
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EXECUTIVE  SUMMARY 


Many  military  communication  systems  are  required  to  operate  in  environments 
where  an  enemy  may  attempt  to  interfere  with  the  corrununication.  One  of  the  most  serious 
threats  is  the  use  of  jammers.  This  threat  has  caused  designers  of  communications  systems 
to  counter  with  the  use  of  spread  spectrum  techniques  such  as  frequency-hopping. 

Frequency-hopped  spread  spectrum  satellite  commimication  systems  provide  robust 
communications  resilient  to  many  types  and  levels  of  jamming.  Unfortunately  this  increase 
in  resilience  is  offset  by  an  increase  in  complexity  while  establishing  the  communication 
link,  termed  synchronization.  In  order  for  the  entire  communication  system  to  be  effective 
against  hostile  and  inadvertent  interference  the  commimication  link  and  synchronization 
process  must  be  protected  to  the  same  level. 

This  document  examines  a  two-step  downlink  synchronization  method.  The  method 
involves  the  transmission  by  the  satellite  of  synchronization  (sync)  hops  that  are  modulated 
with  a  synchronization  code  and  are  received  by  the  ground  terminal  where  the  two-step 
algorithm  is  implemented.  The  first  step  is  the  envelope  detection  of  the  sync  hops,  coarse 
time  synchronization.  When  the  terminal  is  conHdent  that  the  coarse  synchronization 
process  is  complete  the  second  step  is  activated,  fine  time  synchronization.  Fine  time 
synchronization  involves  correlating  a  downlink  sync  hop  with  an  internally  generated 
reference  sync  hop  to  obtain  an  accurate  time  estimate.  Additionally,  a  prediction  of  the 
frequency  offset  can  be  made  to  fine  tune  the  local  oscillator. 

Computer  simulations  were  made  of  the  synchronization  process  using 
synchronization  codes  of  varying  lengths  to  verify  the  algorithm.  These  simulations 
demonstrate  that  the  proposed  algorithm  is  capable  of  operating  under  severe  jamming 
conditions  and  is  compatible  with  the  corresponding  communication  link. 
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1.0 

INTRODUCTION 


Frequency-hopped  (FH)  spread  spectrum  satellite  communication  systems  provide  a 
robust  communications  link  resilient  to  both  hostile  and  inadvertent  interference. 
Unfortunately  the  increase  in  resilient  is  offset  by  an  increase  in  complexity  of  the  initiation 
of  the  communication  link,  i.e.  the  synchronization  of  the  ground  terminal  with  the  satellite. 
If  the  system  is  to  perform  satisfactorily  in  a  hostile  environment  the  synchronization 
process  must  be  at  least  as  robust  as  the  communications  link. 

This  document  describes  the  simulation  of  a  downlink  synchronization  process 
proposed  for  a  FH  communications  system.  An  explanation  of  the  synchronization 
algorithm  is  given  in  the  section  2  followed  by  a  description  of  the  simulation  process  in 
section  3.  The  results  of  the  simulation  process  are  presented  in  section  4  with  a  conclusion 
of  the  results  in  section  5. 


2.0 


SYNCHRONIZATION  PROCESS 


This  section  describes  the  synchronization  process  used  in  this  document  The 
downlink  synchronization  process  involves  spatial,  temporal,  and  frequency 
synchronization.  The  spatial  synchronization  process  is  outside  the  scope  of  this  document 
and  is  assumed  to  have  been  achieved  before  initiating  temporal  and  frequency 
synchronization.  The  equations  used  in  the  determination  of  die  coarse  time  estimate,  Hne 
time  estimate,  and  frequency  offset  are  derived.  In  addition  lower  bounds  for  the  coarse  and 
fine  time  estimate  are  developed. 

The  downlink  synchronization  process  (for  the  system  under  discussion)  is  done 
within  the  ground  terminal  with  the  aid  of  synchronization  (sync)  hops  transmitted  by  the 
satellite.  The  terminal  starts  by  detecting  the  sync  hops  and  calculating  a  coarse  time 
estimate  (coarse  time  synchronization),  followed  by  calculating  a  fine  time  estimate  (tine 
time  synchronization).  When  downlink  time  synchronization  is  complete,  an  estimate  of  the 
frequency  offset  is  made.  The  estimate  of  frequency  offset  can  be  used  by  the  downlink 
receiver  to  tine  tune  the  local  oscillator. 

The  sync  hops  are  modulated  using  phase-shift  keying  (PSK)  with  a  code  chosen 
for  its  good  autocorrelation  properties.  That  is,  the  code  has  the  lowest  maximum  sidelobe 
for  the  given  code  length.  Tt^  different  lengths  were  chosen:  a  short  13-bit  Barker  code,  a 
medium  length  28-bit  code,  and  a  long  1 12-bit  code.  The  28-bit  and  1 12-bit  codes  were 
found  by  using  a  computer  search  technique  to  select  the  best  code.  The  codes  used  in  this 
document  are  given  in  Table  1.  The  autocorrelation  plots  for  these  codes  are  given  in 
Rgures  lA,  IB,  and  1C. 


Code  length 

Code  sequence 
(hexade^al) 

13  bits 

159F 

28  bits 

18FD  5B6 

1 12  bits 

FA74  2D84  1B42  E5F3  8C95  D5D4  98E7 

Table  1  Synchronization  codes  for  various  code  lengths. 

Sync  hops  are  transmitted  by  the  satellite  in  bursts.  A  sync  burst  consists  of  4  sync 
hops,  each  separated  by  one  hop.  The  structure  of  a  sync  burst  is  shown  in  Figure  2. 
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Figure  1 A  Autocorrelation  plot  for  13-bit  code. 
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Figure  2  Structure  of  synchronization  burst. 


Sync  Hop  Detection  Algorithm 

The  first  step  in  achieving  synchronization  is  to  detect  the  sync  burst  (envelope 
detection).  Referring  to  Figure  3,  the  detection  process  is: 

a.  the  terminal  dehopping  synthesizer  is  set  to  the  frequency  of  the  first  hop  of 
the  synchronization  burst 

b.  the  downlink  signal  is  dehopped,  bandpass  filtered,  squared,  and  lowpass 
filtered. 

c.  the  resultant  analog  signal  is  sampled  (3  samples  per  hop)  and  converted  to  a 
digital  value. 

d.  the  digital  sampled  value  is  added  to  the  two  previous  samples  and  compared 
to  a  predetermined  threshold. 

e.  when  the  threshold  is  exceeded  the  synchronization  process  jumps  ahead  and 
looks  for  the  next  sync  hop. 

f.  if  the  threshold  is  exceeded  for  three  consecutive  sync  hops  the  sync  hop  burst 
is  assumed  to  be  detected  and  the  fourth  sync  hop  in  the  series  is  used  to 
calculate  the  coarse  timing  estimate.  If  a  sync  hop  is  not  detected  in  the  second 
or  third  hop,  the  process  resets  (i.e.  dehopping  synthesizer  is  set  to  the 
frequency  of  the  first  hop)  and  continues  the  search. 


Coarse  Time  Estimate  Calculation 

Once  the  sync  burst  is  detected,  the  fourth  hop  of  the  burst  is  used  to  calculate  the 
coarse  time  estimate.  The  output  from  the  detection  algorithm  traces  out  the  autocorrelation 
triangle  of  the  sync  hop  envelope.  The  peak  of  this  triangle  indicates  the  point  of  zero  timing 
offset.  The  base  of  the  triangle  is  two  hop  intervals  wide  (see  Figure  4).  The  timing 
mismatch  from  the  largest  autocorrelation  is  given  by  (1),  which  is  derived  from  the 
geometry  of  the  triangle  [1,  equation  (6)]. 


f  V2tj)op  •  Vi(thoD  -  ts) 
toff-  V1  +  V2 


where 


thop  =  width  of  hop 

tg  =  time  between  two  samples 

V 1  =  amplitude  of  largest  autocorrelation  value 

V2  =  amplitude  of  second  largest  autocorrelation  value 
ti,  t2  =  time  of  sample  V i,  V2 


(1) 


5 


band  pass  squanng  low  pass  a  to  u 

filter  circuit  filler  convertor 


Figure  3  Block  diagram  of  sync  hop  detection  and  coarse  time  estimate. 
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Figure  4  Autocorrelation  triangle  of  the  sync  hop  envelope. 


Coarse  timing  uncertainty  bound 


Employing  techniques  first  used  in  radar  to  determine  the  accuracy  of  radar 
measurements,  a  lower  bound  on  the  accuracy  of  the  timing  estimate  can  be  derived.  This 
bound  is  usually  referred  to  as  the  Cramer-Rao  (CR)  bound.  The  general  time-delay  rms 
error  lower  bound  (5Tr)  is  defined  by  [2,  equation  1 1.17] 


where:  p  =  effective  bandwidth, 

E  =  signal  energy,  and 

No  =  noise  spectral  density. 


To  derive  the  lower-bound  for  coarse  for  our  simulation,  we  first  make  the 
assumption  that  the  signal  is  a  bandwidth-limited  rectangular  pulse  [2,  equation  11.18]  to 
obtain 


_ nBX-sintcBT _ 

Si(jcBT)  +  (costcBT-IVtcBt 


(3) 


where 


Si(nBT)  = 


tcBT 


du,  the  sine  integral  function 


B  5  bandwidth  of  filter 

T  =  pulse  width 


The  normalized  timing  uncertainty  bound  can  de  derived  by  dividing  both  sides  of  (3)  by  T; 
replacing  the  energy  E  by  NEb  (N  is  the  number  of  bits  in  the  hop  and  Eb  is  the  energy  of 
one  bit);  replacing  T,  the  hop  duration,  by  NXb  (Tb  is  the  duration  of  one  bit);  and  re¬ 
arranging  terms  giving 


8Tr  1 


where  P  is  derived  from  (3). 


(4) 


Equation  (S)  is  the  normalized  timing  uncertainty  bound  for  a  continuous  wave  (CW)  pulse. 
For  our  simulation,  the  Eb/No  ratio  should  be  replaced  by  an  effective  Eb^o  which  includes 
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losses  due  to  PSK  modulation.  The  SNR  for  a  PSK  modulated  signal  of  the  lowpass 
component  of  the  signal  is  given  in  [4]  as 


SNRpsK  = 


NppsK^ 

+  2rpsK) 


wheie:  Ppsk  =  Eb^No  for  a  PSK  modulated  pulse 

4  =  2,  the  time  bandwidth  product 

The  SNR  for  an  equivalent  CW  pulse  with  same  ^  can  be  derived  as 

SNRcw  =  (7/4  +  2NPcw) 

where:  P(^  =  Eb/No  for  a  CW  pulse 


(6) 


Using  the  same  technique  as  [4],  an  expression  can  be  derived  that  transforms  the  Et/No  for 
a  PSK  modulated  pulse  to  the  effective  Eb/No  for  a  CW  pulse.  This  transform  is  given  by 


f 


CW  “  2+2ppsK 

Therefore  the  lower  bound  for  normalized  coarse  timing  uncertainty  is: 
STr,  1 

/“t 

jq  Pew 


(7) 


(8) 


Fine  Time  Estimate  Calculation 

After  a  coarse  time  estimate  has  been  obtained,  fine  time  estimation  is  done  on  the 
next  sync  burst  Referring  to  Hgure  5,  the  time  estimation  process  is: 

a.  the  terminal  dehopping  synthesizer  is  set  to  the  frequency  of  the  first  sync 
hop  in  the  next  sync  burst 

b.  the  downlink  signal  is  dehopped  and  bandpass  filtered. 
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c.  the  hop  is  demodulated  by  mixing  the  signal  with  the  expected  offset 
frequency  (f  o).  The  observation  time  starts  1/2  a  hop  before  tte  coarse  time 
estimate  and  extends  1/2  hop  past  the  estimate  stop  time.  This  allows  for  a 
possible  coarse  time  estimate  error  of  ±  1/2  hop. 

d.  the  resultant  signal  is  sampled  (2  samples  per  code  bit). 

e.  the  digital  sample  is  cross-correlated  with  the  sync  code  to  obtain  a  matched 
filter  output  As  with  the  coarse  timing,  a  triangle  is  traced  with  the  peak 
indicating  the  point  of  zero  timing  mismatch  with  the  base  two  bits  wide. 

Calculation  of  the  timing  offset  is  derived  from  the  same  geometry  as  (1)  and  is  given  by 
(2).  The  fine  time  estimate  of  to  is  given  by 


*  _  V2th-Vi(th-t.) 
^  “  Vi  +  V2 


(9) 


where 


tb  =  width  of  a  bit 

ts  =  time  between  two  samples 

V I  =  amplitude  of  largest  crosscorrelation  value 

V2  =  amplitude  of  second  largest  crosscorrelation  value 


Fine  timing  uncertainty  bound 

To  derive  the  fine  timing  estimate  lower  bound,  we  make  use  of  the  CR  bound  and 
calculate  the  value  of  P  from  (3)  using  numerical  integration  of  [2,  equation  1 1.16] 


-00 

J  (2)t/)2  [S(/)]2  if 

.  . ao) 

/  [S(/)]2  if 

00 

where:  S(f)  =  is  the  Fourier  transform  of  the  input  signal 

s(t). 


The  value  of  P  is  dependent  upon  the  bandwidth  and  shape  of  the  input  filter.  The  values  of 
Pt  (fme  time  synchronization)  for  the  filters  and  pulse  used  in  the  simulations  are  given  in 
Table  2.  The  lower  bound  for  nonnalized  fine  timing  uncertainty  is: 


Pulse 

duration 

tb 

(hop) 

Time 

bandwidth 

product 

px 

Filter  comer 
firequency 
(Hz) 

1/13 

1 

1.882 

6.5 

1/13 

2 

2.201 

13.0 

1/13 

4 

2.979 

26.0 

1/13 

6 

3.685 

39.0 

1/28 

2 

2.149 

28.0 

1/28 

4 

2.979 

56.0 

1/28 

6 

4.262 

112.0 

1/56 

2 

2.139 

56.0 

1/56 

4 

2.971 

56.0 

1/112 

2 

2.135 

112.0 

J 


Table  2  Values  of  Px  for  various  time  bandwidth  products  (fine  time  sync). 


Offset  frequency  calculation 

An  estimate  of  ftequency  error  can  be  made  by  comparing  the  matched  filter  outputs 
of  the  sync  hop  at  0,  fs  and  -fs,  where  fs  is  chosen  to  be  30%  of  the  hop  rate  for  the 
simulations  described  in  this  document  The  frequency  estimate  is  approximated  by 


est  = 


(Pfs  -  P-fs)  X  X  /Correction 

Pfs  +  P-fs  +  PO 


whrae 


Pk  =  |Xk 


(12) 


» 


.  OO 

Xk=  lx(t) 

•OO 

x(t)  =  input  signal  with  coding  removed 

/coirection  =  experimentally  derived  value  to  minimize  estimate  bias 
k  =  -fs.0,orfs 


The  normalized  lower  CR  bound  for  the  frequency  estimation  method  is  obtained 
from  [1,  equation  16]  as 


S/o  0.433 
Rh 
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3.0 


COMPUTER  SIMULATION 


The  synchronization  process  was  simulated  using  a  Monte  Carlo  method  written  in 
FORTRAN.  The  simulations  were  run  on  a  VAX  computer  using  the  VMS  operating 
system  and  on  a  Macintosh.  The  simulation  software  was  divided  into  3  main  segments; 
transmitter,  channel,  and  receiver.  Each  of  these  segments  is  described  in  the  following 
paragraphs. 


Transnnitter 

The  transmit  segment  simulates  the  function  of  the  payload  in  forming  the  sync 
hops  and  transmitting  them  to  the  user  terminals.  The  sync  hop  waveform  is  modulated  with 
a  sync  code  using  binary  PSK  (BPSK)  signal  modulation.  The  sync  code  is  one  of  the 
thiw  codes  described  in  Table  1.  The  BPSK  signal  is  then  upconverted  by  mixing  with  the 
carrier  frequency  plus  some  optional  frequency  offset  The  frequency  offset  represents  a 
mix  of  oscillator  offset  and  Doppler  frequency  offset  calculation  errors.  The  offset  is  a  user 
selected  parameter.  The  frequency  hopping  portion  of  the  system  is  not  simulated;  however 
the  initial  phase  of  each  "hop"  is  a  uniform  random  variable  from  0  to  Zk  . 

The  sync  burst  is  a  digital  representation  of  Figure  2  and  is  7  hops  in  length.  The 
computer  representation  is  a  two  dimensional  array  with  one  dimension  being  the  amplitude 
of  the  digital  signal  and  the  other  dimension  being  sample  number.  This  sync  burst  is 
placed  into  a  larger  array  with  a  random  starting  sample  number.  The  random  start  time  is 
the  actual  timing  offset,  to.  This  larger  array  represents  the  total  time  for  a  synchronization 
attempt  (a  synchronization  attempt  is  defined  in  section  4). 


Channel 

The  channel  segment  simulates  the  waveform  being  perturbed  by  suming  additive 
white  Gaussian  noise  (AWGN)  to  the  output  array  from  the  transmit  section  of  the 
program.  The  power  of  the  signal  is  fixed  while  the  power  of  the  AWGN  is  variable.  The 
ratio  of  the  two  determines  the  Et/No  for  the  simulation.  The  AWGN  is  generated  from 
normally  distributed  random  numbers  using  the  GASDEV  function  from  [3]. 
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Reodver 


The  receive  segment  simulates  the  downlink  synchronization  function  of  the  ground 
terminal.  The  sync  hop  plus  noise  is  first  downconveited  and  then  split  into  the  coarse  and 
fine  synchronization  branches.  Each  signal  is  then  bandpass  filtered  with  a  user  specified 
infinite  impulse  response  (IIR)  filter. 

The  thresholds  used  for  coarse  synchronization  on  the  first  3  hops  are  user 
selectable  and  for  this  report  are  identical.  The  number  of  coarse  samples  per  hop,  the 
number  of  fine  samples  per  bit,  and  the  number  of  crosscorrelations  values  used  in  die  fine 
time  estimate  are  also  user  selectable.  The  values  for  the  simulation  results  described  in 
section  4  are: 

Number  of  coarse  samples  per  hop  3 

Number  of  fine  samples  per  hop  2 

Number  of  matched  filter  values  per  estimate  1 


To  increase  the  accuracy  of  the  estimate  of  the  timing  offset,  to,  the  number  of  fme 
or  coarse  samples  per  bit  could  be  increased.  This  would  have  the  effect  of  more  accurately 
tracing  the  crosscorrelation  peak. 
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4.0 


SIMULATION  RESULTS 


The  results  of  the  computer  simulations  of  the  synchronization  process  are 
described  in  this  section.  Each  computer  run  consisted  of  a  user  specified  number  of 
synchronization  attempts.  An  attempt  is  classified  as  successful  (sync  burst  found),  false 
alarm  (software  not  synchronized  to  a  sync  burst),  or  unsuccessful  (synchronization  was 
not  obtained  in  10  sync  hop  bursts).  As  expected  the  number  of  false  alarms  and 
unsuccessful  attempts  increased  as  the  Ei/No  decreased.  To  insure  the  statistical  validity  of 
each  computer  run  a  minimum  of  500  successful  synchronizations  were  collected. 

Coarse  Hme  Estimates 

Simulations  to  determine  coarse  time  estimates  were  done  for  three  code  lengths: 
13  bits,  28  bits,  and  1 12  bits.  The  RMS  timing  uncertainty  (in  bits)  from  the  simulations  are 
compared  to  the  CR  bound ,  refer  to  (7).  The  results  ate  shown  in  Figures  6, 7,  and  8  versus 
Eb/No.  The  agreement  between  the  simulation  and  bound  for  high  Eb/No  is  very  good. 
However  for  lower  values  of  Eb/No  the  agreonent  deteriates.  This  is  to  be  expected  because 
of  the  assumption  of  large  Eb/No  in  the  derivation  of  the  lower  bound. 

One  important  result  from  the  simulation  is  that  the  synchronization  algorithm  was 
able  to  function  at  low  values  of  Eb/No.  In  all  cases  the  algorithm  was  able  function  down  to 
an  Eb/NoOf  IS  dB.  However,  at  times  the  algorithm  had  to  cycle  many  times  before  a  coarse 
estimate  was  obtained. 

Fine  Hme  Estimates 

Simulations  to  determine  fme  time  estimates  were  done  for  the  same  three  code 
lengths  (13  bits,  28  bits,  and  1 12  bits)  as  for  coarse  time  simulations.  A  filter  with  a  time- 
bandwidth  product  of  2  was  used.  The  RMS  timing  uncertainty  (in  bits)  from  the 
simulations  compared  to  the  CR  bound  are  given  in  Figures  9,  10,  and  1 1  for  estimates 
based  on  one  hop.  The  RMS  error  for  timing  uncertainty  estimates  based  on  a  four  hop 
average  ate  shown  in  Hgutes  12, 13,  and  14. 
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Figure  10  Fine  timing  accuracy  for  a  28-bit  synchronization  code  using  one  hop. 


Figure  1 1  Fine  timing  accuracy  for  a  1 12-bit  synchronization  code  using  one  hop. 
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Figure  12  Fine  timing  accuracy  for  a  13-bit  synchronization  code  using  4  hops. 
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Figure  13  Fine  timing  accuracy  for  a  28-bit  synchronization  code  using  4  hops 
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Figure  14  Fine  timing  accuracy  for  a  1 12-bit  synduonization  code  using  4  hc^. 
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As  in  the  simulation  of  coarse  time  estimate,  the  fine  time  estimation  algorithm  was 
shown  to  function  at  low  values  of  Eb/No.  In  fact  the  fine  algorithm  worked  to  lower  values 
of  Eb/No  than  the  coarse  algorithm.  Figures  9  and  12  (code  length  =13)  show  the  upper 
limit  of  where  the  synchronization  error  is  considered  to  be  acceptable,  less  than  1/20  *  of  a 
bit.  Using  a  single  hop  for  the  fine  time  estimate,  the  algorithm  showed  that  acceptable 
timing  estimates  could  be  achieved  down  to  an  Eb/No  of  8  dB.  If  a  4-hop  average  is  used, 
the  algorithm  was  acceptable  down  to  an  Eb/No  of  4  dB. 

The  improvement  of  the  4-hop  estimate  over  the  1-hop  estimate  approaches  the 
theoretical  3  dB  at  low  E/No,  but  reduces  to  0  dB  as  the  Eb/No  increases  (see  Figure  15). 
The  error  in  timing  estimates  is  mainly  due  to  two  factors;  the  error  caused  by  the 
approximation  of  the  cross  correlation  function  as  a  triangle  and  the  error  caused  by 
variations  of  the  cross  correlation  value  due  to  noise.  At  high  Eb/No  the  main  error  source  is 
due  to  the  approximation  of  the  correlation  function,  which  is  the  same  for  the  1-hop  and  4- 
hop  estimate.  Therefore  no  improvement  is  possible  by  averaging.  At  low  Eb/No  the  main 
source  of  error  is  due  to  noise  and  the  theoretical  3  dB  improvement  can  be  realized  by 
using  the  4-hop  average  (assuming  the  burst-to-burst  timing  offset  is  random).  It  should  be 
noted  that  the  error  due  to  a  poor  approximation  of  the  autocorrelation  function  is  less  than 
0.02tb  for  all  cases  investigated. 

Frequency  O^et  Estimates 

Simulations  to  determine  the  accuracy  of  the  frequency  offset  calculations  were 
done  for  the  same  code  lengths  as  above  (13  bits,  28  bits  and  1 12  bits).  Frequency  errors  of 
3%,  6%,  9%,  and  12%  of  the  hop  rate  were  introduced  into  the  simulation.  Figure  16  shows 
the  resultant  RMS  error  between  the  actual  frequency  offset  and  the  calculated  offset  for  a 
code  length  of  13  bits.  The  RMS  error  is  expressed  as  a  ratio  of  the  hop  rate.  Figure  16 
shows  the  RMS  error  using  only  1  hop  in  the  calculation.  Figure  17  shows  the  RMS  error 
using  a  4-hop  average  in  the  calculation.  For  Eb/No  values  greater  then  5dB  the  errors  for 
the  various  frequency  offsets  are  closely  grouped  around  the  average. 

The  average  RMS  errors  versus  Eb/No  are  given  in  Figure  18  for  all  3  code  lengths. 
These  errors  are  compared  to  the  CR  bound  derived  from  equation  11.  Figure  18  shows  that 
the  RMS  error  for  the  4-hop  average  is  3  dB^  better  than  using  the  1-hop  value  for  all  cases 
as  would  be  expected.  For  RMS  errors  of  less  then  3%  the  bit  error  rate  (BER)  is  degraded 
by  a  negligible  amount  Using  this  criteria  the  estimate  is  sufficiently  accurate  down  to  an 
Efa/No  of  24  dB  using  the  1-hop  value  and  accurate  down  to  an  Efa/No  of  18  dB  the  4-hop 
value.  Table  4  shows  the  accuracy  cutoff  value  for  the  various  code  lengths  in  terms  of 
Eb/No.  If  necessary  the  estimate  could  be  improved  by  using  more  hops  in  the  frequency 
offset  estimate. 


^The  increase  of  accuracy  at  high  Eb/No  is  defined  to  bel/Vn,  where  for  our  case  n  =  4. 

^Analysis  of  an  imperfect  carrier  reference  for  a  BPSK  signal  in  [5]  shows  that  a  3%  carrier 
frequency  error  has  a  minimal  effect  on  the  bit  error  rate. 


19 


Ratio  (dB) 


Figure  IS  Ratio  of  4-hop  fine  timing  estimate  over  1-hop  flue  timing  estimate. 
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ngure  16  Frequency  offset  calculation  accuracy  using  one  hop. 
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Figure  17  Frequency  offset  calculation  accuracy  using  4  hops. 
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Eh/No 

_ EbffJo  1 

13  Bit  Code 

28  Bit  Code 

112  Bit  Code 

1  Hop 

Estimate 

24dB 

13dB 

lOdB 

4dB 

4-Hop 

Estimate 

18dB 

7dB 

4dB 

-2dB 

Table  4  Accuracy  cutoff  value  for  3  code  lengths. 


The  effects  of  temperal  and  frequency  errors  are  graphically  depicted  in  the 
ambiguity  diagrams  for  the  three  synchronization  codes  (Figures  19,  20,  and  21)  at  high 
Eb/No.  As  can  be  seen  in  the  figures  the  sidelobes  increase  as  frequency  offset  increases. 
As  these  sidelobes  increase,  the  probability  that  the  synchronization  algorithm  will  choose 
the  wrong  offset  also  increases.  Comparing  the  three  figures,  it  is  apparent  that  the  longer 
the  synchronization  code,  the  better  its  performance  is  against  frequency  offset  errors. 
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Figure  19.  Ambiguity  Diagram  for  13-Bit  Sync  Code 
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Figure  20.  Ambiguity  Diagram  for  28-Bit  Sync  Code 


5.0 

CONCLUSION 


This  document  has  described  an  method  for  downlink  synchronization  that  is  at 
least  as  robust  as  the  spread  spectrum  communications  system.  The  associated  algorithm 
was  used  with  synchronization  codes  of  three  different  lengths:  a  short  code  (13  bits),  a 
medium  length  code  (28  bits),  and  a  long  code  (1 12  bits).  The  selection  of  which  code  to  be 
used  would  be  dependant  upon  the  required  accuracy  of  the  downlink  time  and  frequency 
and  the  data  link  waveform  structure.  In  general,  the  longer  the  code,  the  greater  the 
accuracy. 
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